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ABSTRACT 


We list and discuss the tasks we have accomplished and considered. In some areas 
more has been accomplished than others. An extra task was a calculation comparing 
electron mobilities in Mercury Manganese Telluride with Mercury Cadmium Telluride 
given in 1 H. We then list the reports and papers produced and follow that with either 
abstracts or the papers themselves. In one key paper we obtain good results between 
experiment and theory in Mercury Zinc Telluride and also find it typically has mobilities 
competitive with Mercury Cadmium Telluride. In the appendix we have a relatively 
complete set of references. 
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I. LIST OF TASKS ACCOMPLISHED/CONSIDERED 

A. Characterization of Donors and Acceptors 

The study of donors and acceptors in narrow gap semiconductors is by no means 
simple. Misplaced Hg in MCT seems to be a major cause of intrinsic doping. Hg 
vacancies can produce a p type material and extra Hg interstitials an n type. Indium 
doping has caused p type MCT to become n type as expected. A little more generally, 
group I and II elements should be acceptors and donors on Hg, Cd (group II) metal 
sites while V and VII elements should be acceptors and donors on Te (group VI) sites. 
The situation is very complicated, however, and experiment is needed to verify these 
suppositions. Furthermore, the situation is complicated by the fact that we can have 
both deep and shallow states. Vacancies, interstitials and even substitutional 
impurities can cause deep levels. In any case deep levels produced by a short range 
potential, are often not well characterized, and usually produce a level near the middle 
of the gap. If we neglect deep levels, we are still left with shallow levels. Shallow 
states can, at least for small x, merge with the conduction band. We can understand 
this merger as related to the Mott transition. When the density of defects increases, 
discrete levels may become broadened and produce energy states that may join with, 
for example, the conduction band in the case of donors. In fact, it has been predicted 
that Hg.,_ x Cd x Te with a (shallow) donor concentration of at least 10 14 cm* 3 would not 
have bound states and this seems to be well confirmed experimentally. In general 
there appears to be no freezeout of carriers on shallow impurities for bandgaps in, as 
grown, MCT which has x < .2. There are in fact other reasons for the merging with the 
conduction band. Donor levels can be screened by electrons in the conduction band 
and this screening can decrease the binding energy. We have examined this in a 
paper referenced below. In addition, when the binding energy gets small enough the 
donor wave functions may hybridize with the band functions and again we get 
merging. 

Note about shallow donors . 

The simplest theory of shallow impurities characterizes them as hydrogenic. For non- 
degenerate bands with parabolic dispersion relations, the “envelope" function of the 
impurities satisfies; 

+ 2 ? 2 

[ - -2— V - — ] F = E F 
* Kr 


where m* is the effective mass and K is the (relative) dielectric constant. 
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13,600 m 
tB = 2 7 ^ 

For a Is ground state , the binding energy is K ‘ 


mev 


with wave function 


Fls“ 




where 

m 0 i 

a B =0.53xKx — °A 
m 

Our situation has many complications beyond this (Bastard p. 119-120). We don’t 
have parabolic bands and the donors are often not well characterized. Another 
complication arises from the fact that in the c o mpound semiconductors that we 
consider, the lattice is charged so electrons interact with it even at absolute zero. Thus 
the stationary states are polarons. The effect of the lattice can only be neglected for 
tightly bound electrons - not the usual ca se fo r shallow donors. For weak binding one 
must view the situation as a polaron bound in a hydrogenic orbit 

SUMMARY - SHALLOW LEVEL DEFECTS 

These have a long range potential and their wave function spreads over many 
neighboring atoms. 

Have a hydrogen like spectrum as modified by dielectric constant and effective mass. 
Usually have ionization energies much less than the gap energy. 

Connected to the adjacent band (that is Eigenfunction mix with nearest band states). 
Due to screening and hybridization may merge with nearest band. 

They control conductivity. 

Note about deeper levels . 

The nature of donors in Ga-|_ x Al x As can be quite different depending on the x value 
(Weisbuch and Vinter, 47-50). For Si donors the energy levels can change from 6 
meV to 160 meV as one changes from direct to indirect materials. The transition takes 
place at about x = 0.24. 
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Erbarut has treated the electronic structure of ideal vacancies in MCT. His work is 
based on an extension of the Slater-Koster model and does not include lattice 
relaxation. The energy levels obtained are often in the range of deep bound levels. 

At any rate, deep levels are often treated along Slater-Koster lines. Effective mass 
techniques are often not appropriate for deep levels. Hjalmarson has defined a deep 
impurity to be one “whose short-range central-cell potential alone is sufficiently strong 
to bind a state”. 

SUMMARY - DEEP LEVEL DEFECTS 

Short range potential, wave function relatively well localized. 

Not a hydrogen like spectrum - use tight binding approximation. 

Often have ionization energy comparable to half the gap energy. 

Connected to the conduction band and valence band. 

Often hard to even identify a model for center. 


They control lifetime of carriers. 

General note . 

Capper has given a review of impurity behavior in bulk MCT. He points out that 
“extensive doping using elements which provide stability of electrical properties will 
become increasingly important in future infrared devices based on MCT.” A common 
way to study defects is by photoluminescence spectra as discovered by Werner et al. 
Native defects are highly sensitive to deviations in stochiometry and to stress. Point 
defects come in a variety of forms; lattice defects (vacancy, self interstitial and antisite) 
and foreign atom (substitutional and interstitial). 

Donors, acceptors and other defects can be important in lowering the mobility at low 
temperature, but other scattering mechanisms are much more important at higher 
temperatures. Defects determine to a large extent the quality of these materials for 
infrared detectors. 


3a. 
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B. Electron-Longitudinal Optic Mode Scattering 

We have studied the Random Element Isodisplacement Model of two-mode behavior. 
We have also reviewed models for establishing criteria for two mode versus one mode 
behavior. We have studied effective charges and their meaning for scattering 
electrons. We have also reviewed the literature. All of these things are discussed in 
more detail below. 

1. Szigeti, Transverse and Callen Effective Charge 

Ionic polarizability is caused by the displacements of the ions as a whole while 
electronic polarizability is caused by the displacement of the electron cloud 
relative to the nuclei. The displacement of the ions changes the interatomic 
distance which in turn affects the distribution of the electron cloud. This effect is 
a short range interaction, since it only depends on the displacement of the near 
neighbors. There is another type of interaction between atomic and electronic 
polarizations through the internal field. This effect is a long range effect since 
the distant dipoles as well as the near dipoles contribute. 

Szigeti defined a correction factor s for deviation from heteropolar behavior as a 
result of mutual distortions of neighboring ions due to their overlapping. Szigeti 
charge represents the short range interaction of the electronic and atomic 
displacements. One might expect to find s very different from one if the ions 
penetrate considerably into each other or if the bond has partial homopolar 
character and its polarity changes rapidly with the distance between the atoms. 
Szigeti has shown how the factor s can be obtained in terms of values that can 
all be experimentally measured. 

The transverse charge is defined to give the local polarization (electronic and 
atomic) caused by relative displacements. The polarization includes that due to 
ionic motion and that induced in the electronic charges in the absence of 
electric fields. The resulting field which scatters electrons is accounted for by 
separately dividing by e 0 (giving rise to the definition of the Callen charge). 
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The Szigeti effective charge e s = se is used in computing the force on a ion. Distortion 
of the charge affects the value of e s and causes s * 1 . The distortion is due to overlap. 
By Born and Huang 

P = — { sze (u + “ U. ) + (a + + a. ) P 

3 

Assuming E = 0, a ± are the polarizability of the ions and u + - u_ is the relative 
displacement of the ions, z is the charge of the ions and v a is the volume per 

molecule. Solving for p gives us 


P = 


1 sze (u + - u. ) 


v 


a 1 


4n 

3v„ 


(a + + a. ) 


Then using the Clausius-Mosotti equation 


An / s e<» 1 

— — (a + + a. ) = V 

3 eoo + 2 


This means that the polarization can be written 

P = — Uoo + 2) (u + - u. ) 

v a 3 

The quantity 


se = 


is called the transverse charge. Thus 

P = — e T (u + - u_ ) 

V a 

This P produces a field which scatters electrons. But in calculating this field we must 
divide by since it is long range and the field is affected by other electrons. So 
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another relevant charge, called the Callen charge is 

e c = c j /e oo 


SUMMARY 


I . 


e s : The Szigeti Charge e s = se 

This is the effective charge used in figuring the force on an ion. Its value is 
affected by distortion due to overlap. It is a measure of deviation from ideal 
heteropolar behavior. 


e ro + 2 

e T = e s 

er. The Transverse Effective Charge 3 

This is defined so as to give local polarization induced by relative displacement. 
We are interested in induced polarization in the absence of electric fields. 


e c : The Callen Effective Charge 


e c = c j/eoo 


€<» + 2 
3e<n 


Note for rigid ions = 1 and e c = e T = e s . 
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2. One and Two Mode Behavior 

In binary alloy systems, lattice vibration spectra can show two different 
behaviors. The first is “one mode behavior” in which the phonon frequency 
changes continuously from the value that corresponds to one end member to 
that of the other end member with an approximately constant strength. The 
second is “two mode behavior” in which two frequencies occur with a strength 
approximately equal to the fractional composition of each component. 

Experimentally, this behavior can be seen through infrared and Raman 
spectroscopy. A Kramers-Kronig dispersion analysis can be applied to 
reflection spectra to get the real part of the dielectric constant Re e(cj), the 
imaginary part of the dielectric constant Im e(uj), and the imaginary part of the 
reciprocal of the dielectric constant -Im (l/e(oo)). From these functions the 
frequencies of the transverse optical phonons (TO) and those of longitudinal 
optical (LO) phonons are determined . The frequencies of the TO phonons are 
determined from the positions of the maxima of the Im e(co) curve and the 
frequencies of the LO phonons from the positions of maxima of -Im (l/e(oo)) or 
from the points where the Re e(oj) curve crosses the oj - axis. In the case of 
the two mode behavior, the Kramer-Kronig analysis yields for each x - value two 
frequencies of the transverse optical phonons (TO! and T0 2 ) and two 
frequencies of the longitudinal optical phonons (LO! and L0 2 ). 

Now returning to the ideas of one and two mode behavior. Consider the 
pseudobinary compound AB!_ X C X composed of atoms A, B, and C. Two mode 

behavior in AB!_ X C X may also be seen experimentally in reflectivity spectra. 
Supposedly, this means the interactions of a cation with one anion and a cation 
with another anion are independent, thus the AB and AC vibrations are 
independent. Actually one may see one or two mode behavior depending on 
the situation. 

In one mode behavior the frequencies of long wave length phonons vary 
continuously from x = 0 (AB) to x = 1 (AC). 

In two mode behavior one has two sets of phonon frequencies for each x * 0 or 
1 . The set that starts from x = 0 starts from AB type behavior and the set from x = 
1 starts from BC type behavior. The mole fractions of AB and AC determine by 
proportionality the strength of the frequencies. 
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A sketch of the two kinds of behavior looks like: 



The following convention is used in naming modes. An impurity mode lying 
above the optic band of the host lattice is a local mode. Note this occurs at 
x = 1 . A gap mode lies below the band of the host lattice. Note this occurs at 
x = 0. There are also more complicated situations. Obviously one can reverse 
directions by letting x -* x' = 1-x , 1 - x -* x'. 

The actual criteria for one - and two- mode behavior can be somewhat 
complicated. 

Let 


€ = 1 - M'/M 


M' 

M' 


M 

M 


where 
and for 


< 


€ < 0 
0 < e < 1- 
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Insight into the times when we should expect two mode behavior and one mode 
behavior can be gained by considering the vibrational behavior of a one 
dimensional linear diatomic chain. 

If we substitute for the lighter mass with an even lighter impurity (so e > 0) a 
localized mode will rise out of the top of the optical band. If we substitute with a 
heavier mass (e < 0) an optical gap mode falls out of the bottom of the same 
optical branch. If we substitute for the heavier mass with e > 0, then a local 
mode rises from the top of the optical branch and a gap mode rises out of the 
top of the acoustical branch. If one substitutes for the heavier mass with e < 0, 
no new gap or local modes appears^ Obviously, one needs local and gap 
modes for two mode behavior. 

Chang and Mitra give the following conditions: 

I. M b < (j ac so M b < M a and M c (assuming M B < M c for AB.,_ X C x ) 
or 

II. M c <^ ab soM c <M A andM B (assuming M c < M B for AB.). X C X ) 

If (I.) or (II.) is satisfied, then AB-|_ X C X is two mode, otherwise it is one mode. 
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3. The random element isodisplacement model 

There are several models for calculations in mixed crystal systems. In a virtual 
crystal model all masses and spring constants are taken to have the same 
average value. 

Verleur and Barker have accounted for two mode behavior in some mixed 
crystals using a model based on short range clustering. 

This is a good working model proposed by Chen and developed by Chang and 
Mitra to account for two mode behavior in crystals of the form AB.|_ X C X where x 
gives the concentration of C relative to A, See Gorska and Nazarewicz. There 
are two basic assumptions 

(a) All atoms of the same kind move together (isodisplacement). 

(b) We can think of the crystal as having two sublattices. One is a lattice of 
atoms of the type A and the other is a random distribution of atoms of the 
type B and C with relative co ncentrations 1-x and x. 


Let m A , m B , m c be the masses of the atoms, u A , u B> u c be the displacements 
from equilibrium, e A , e B , e c be the effective charges of the ions, E ef{ be the 
effective field, and F AB , F AC , and F BC be the force constants. 

Thus we can write: 


M a U a =-(1 -x) F AB (U A -U B )-xF AC (U A -Uc) + e A E e{f (i a ) 

MbUb = ‘ f ab(Ub'U a )-xF bc (U b -Uc) - e B E efj (ib) 

McUc= ‘ Fac(Uc Ua) - (1 - x) F bc (U c - Us) - e c E eff (i c ) 

We assume a linear dependence of the force constants on the composition: 

Fab = F ac _ F bc = j _ q x 

Fabo F aco F bco (2) 


The effective charges must be determined by the electrical neutrality condition: 
e A -(1 -x)e B -xe c = 0. ( 3 \ 
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We start out by looking at long-wavelength transverse optical phonons for 
which E « 0. For this case 


P = N [ e A U A - (1 - x) e B U B - x e c U c ] + N a E ef{ 


(4a) 


where 


E ff = E + — P = — P 
eff 3 3 


(4b) 


The electronic polarizability is determined by the Clausius-Mossotti equation 


i2. N« - £ “ ' 1 


* 2 

Using Eqs. (4a), (4b) and (5) we find: 

Coo + 2 


(5) 


P = 


N [ e A U A - (1 - x) e B U R - x e c U c ] 


B u B 


( 6 ) 


and 


Eeff - 42- N [ e A U A - (1 - x) e B U B - x e c U c ] 

3 3 (7 ) 


Using Eq. (3), Eq. (7) becomes: 


E e „ - -y- N [ (1 - x) e B (U A - U E ) * x e c (U A - U c ) ] 

Repeated use of Eq. (3) allows us to express E eff as involving (U A - U B ) , 
(U B -U C ), or (U C -U A ),(U C -U B ). 


M A U A = - (1 - x) F’ ab (U A - U B ) - x F' ac (U a - U c ) 
M b Ub - * F' ab (U a - U B ) - x F' bc (U b - U c ) 

M c iic - F' ac (U a - U c ) * (1 - x) F' bc (U b - U c ) 


(8a) 

(8b) 

(8c) 


13 . 


where 


F' 


AB " r AB 


F'ac= Fac- 

F'bc = Fbc + 

Now defining 

Wi=U A -U B 
W 2 = U B - U c 
W 1 + W 2 = U A -U C 
The equations become: 

W [ = ~ ( K i + K j 2 ) W j _ K ] 2 V/ 2 


4n 

€00 + 2 

N e A e B 


3 

3 

(9a) 

4rc 

£co + 2 

N e A e c 


3 

3 

(9b) 

4rt 

€„ + 2 

N e B e c 


3 

3 

(9c) 


W 2 = (K, - K 2 + K , 2 - K 21 ) W, - (K 2 - K 12 ) W 2 


(10a) 

(10b) 

(10c) 


(lla) 

(llb) 


where 


K 


= (1 - x) Iab + Iab + x F'bc 


M 


M 


B 


M 


B 


K 2 = x -^-2 + Hm + (1 - x) — ^ 
M a M c M c 


(12a) 


(12b) 


K,2 - x 


F' ac _ „ F’bc + q _ x ) F'bc 


M 


M, 


Mi 


(12c) 
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K 2| = (1 - x) 1*1 - (1 . X ) £bc 
m A m c 


(1 2d) 


Thus 


oj 4 - CK ! + K 2 ) w 2 + (KjK 2 - K , 2 K 2 1) = 0 


(13) 


We now redo this for the case of longitudinal optical phonons with e = 0, or the 
more familiar form E = - 4 ttP. This is given by Kittel (p. 275) and follows from 


P = — — - E with e = 0. Using 

4n 


E = - 4kP 


(14) 


E„f,= E+^P 


— P 
3 ■ 


(15) 


Combining Eqs. (15), (5) and (4a), we find: 


P = — - (e A U A " (1 - x) e B U B - x e c U c ) 


(16a) 


Eeff = - € °1 + 5 (e A U A - (1 - x) e B U B - x e c U c ) 

3 3e„ 


(16b) 


Thus we find 


M A U A = - (1 - x) F ab (U a - U B ) - x F ac (U a - U c ) 


M b Ub = -F ab (U b -U a )-xF bc (U b -U c ) 


(17a) 

(17b) 


McU c =-F ac (U c -U a )-(1-x) F Rn (U c - U B ) 


BC 


(17c) 
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where 


F ab = F ab * M N e A e B 

(18a) 

f ac = F*c*yN e A e c 

^ 0€ oo 

(18b) 

tT it 8 TC t /^oo + 2. 

F bc ■ F bc - 3 N e B e c ( ^ ). 

(18c) 

Thus we have the same equation for transverse and longitudinal optic phonons 

except that F AB replaces F AB . 


We assume that N and e M depend on x as follows: 


N = (1 - x) N ab + x N ac 

(19a) 

G 00 ~ 0 ^00 ) AB ^ G 00 ) AC 

(19b) 

and we also assume 




(20a) 


e c = e AC (20b) 

★ ★ 

so e A = (1 - x) e AB + x e AC (20c) 

where e* is the Szigeti effective charge. Now it is a matter of solving Eq. (13). 

Summary for crystals of the form AB 1 . X C X , we need to solve Eq. (1 3) for oj 2 with 
the K.’s defined by Eq. (12) in which the F' ’ S are replaced by F" ’ S with Eq. 

(20) used. The solution gives oj 2 as a function of x. 

For the spherical case of long wavelength optical phonons at x = 0 and x = 1 , 
we find the following initial conditions. 
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X = O 


2 _ F ABO A 8n ^AB ( e AB) 

W L 0 AB “7“ 


( e °° , A A + ^ 


FAB 


fab 


3€ 


00 , AB 


2 _ f aco 

^ga p, AB 


+ F 


BCO 


M ( 


X = 1 


2 _ F aco 0 ~ ®) 8n FI A c (e AC) 

w lo ac + -=- ~ 


e °° , A C 


+ 2 \ 


FAC 


F A C 


\ 00 } A C / 


2 F AB0 + f bco .j _ ^ 


^local, AC 


M 


B 


where p AB M A + M B an p AC M A + M C ' 

There is also the question as to whether the two mode theory for photon 
scattering is the same as the two mode theory needed for electron scattering. 
One might say that the Green’s function for scattering in the two cases would be 
clearly different. However we are interested in long wavelength modes mostly, 
and at any rate the Green’s function probably depends more on frequency than 
wave vector. Thus we assume our approximation is not a bad one. This is 
especially true because we assume each optic mode only has one frequency 
which is independent of wave vector. 
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In the following graph we plot MCT mobility versus temperature for several different 
values of transverse charge which have appeared in the literature. The main 
difference appears at intermediate temperature where variations of 10% or so can be 
seen. See also the accompanying table where numbers and references to the values 
of transverse charge are given. It seems clear that the effective charge determining 
the strength of scattering is not rigorously determined. However, different reasonable 
choices seem to have little effect on the mobility vs. temperature. 


Electron mobility (cm2 /Vj) 
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Effect of using different values of 
transverse charge on electron 


mobility 
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C. Review of Alloy Disorder Scattering 

The alloy Hamiltonian can be written, in an approximation, as a uniform virtual 
Hamiltonian with a residual alloy disorder potential. The virtual crystal approximation 
is used in alloys which are formed by ternary solid solutions between ll-VI binary 
compounds. In the virtual crystal approximation (VCA) the actual non periodic 
structure is replaced by an average which is periodic. We consider a crystal of the 
form A X B-,_ X C. The lattice of the atoms A is fee and the lattice of the B, C atoms 
randomly occupy a second fee lattice. Thus in VCA, the random potential, created by 
B(V b ) and C(V c ) atoms is not used. Instead, a periodic potential which is the weighted 
averages of these two is appropriate. In equation form, we write 


< V > = V A + x V c + (f - x) V B . 


The use of the average potential leads again to Bloch waves and effective masses, but 
the residual disorder then leads to scattering between the Bloch waves and to a 
compositional dependence of the band gaps. Only when the two binaries making up 
the pseudobinary alloy have very similar band structure and lattices, does the band 
gap vary approximately linearly with composition. After analysis, this typically leads to 
a scattering term proportional to x(1-x), however things are really not so simple as one 
needs to know with what strength to treat this scattering term. The answer that is 
usually given is that the strength is proportional to the difference between two band 
widths. For example in MCT, the difference in the bandwidth between HgTe and CdTe. 
Kossut has looked at the scattering in more detail and finds that it actually depends on 
three matrix elements involving the difference in potential between Hg and Cd. The 
one involving the spin-orbit interaction is apparently zero for MCT, the other two, are 
not necessarily zero. Only when two of th e th ree are zero is the scattering proportional 
to the difference in the band gap energies from HgTe and CdTe. One further problem 
is the non uniformity of the disorder via, for example, clustering. 
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D. Other Scattering Mechanisms 

Obviously dislocations are imperfections in the crystal and as such they can act as 
scattering centers. For example, edge dislocations create strain fields which can 
scatter electrons due to coupling to the deformation potential. Dislocations can also 
become charged and thus act as scattering centers due to their charge which may be 
screened by free carriers or impurities or both. A discussion of dislocation scattering 
has been given by Zawadzki. It should be pointed out that this subject is not only of 
academic interest. Pelliciari, Destefanis, and DiCioccio have found evidence of 
anomalous behavior in low n type mercury cadmium telluride induced by extended 
defects. Their samples have shown a kink in the mobility curve between 77 and 100 
K. Transmission Electron Microscopy experiments have correlated this kink with the 
appearance of dislocations and a significant decrease in the mobility. A very serious 
practical problem is the effect of temperature gradients on dislocation formation during 
crystal growth. 
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E. Further Work in MZT 

The best summary of further work on MZT that we have done has been summarized in 
a recent paper by us which appeared in the August issue of the Journal of Materials 
Research. One of the main results of t his p aper is a comparison of our work with 
experimental work done by others. We should also mention that although MZT has 
several advantages as a IR detector (over MCT) primarily because of its increased Hg 

stability it also suffers from a greater difficulty in growing uniform crystals. This ’ s 

dEg 

because less Zn than Cd is needed for a given E and dx for E’s of interest j S 

9 

greater for MZT. Thus small errors in x are magnified in energy gap fluctuations. 
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F. Experimental Data 

Experimental data on narrow gap materials other than MCT is relatively sparce. 
i Besides the references given in E we note the papers by M. Voltz et al noted in B. 


We have also “invented” some materials by c hanging the parameters of MZT and MCT 
by factors which are arbitrary but not totally unreasonable. The idea was to check if a 
more or less random but reasonable set of parameters would fit experiment as well as 
any. Happily they did not for otherwise we would have to say that our good results 
were fortuitous. 
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G. Superlattices 

We have several comments to make about superlattices. Unfortunately experimental 
studies require sophisticated equipment such as a molecular beam epitaxy apparatus. 

A superlattice (SL) which has alternating layers of CdTe and HgTe has been proposed 
by Schulman and McGill as a new IR material which can be made by molecular beam 
epitaxy. The first successful molecular beam epitaxy growth of a Hg-Te superlatttice 
occurred in 1 982. The band gap depends on the width of the quantum wells as well 
as the barrier width and height. Thus the formation of HgTe/CdTe superlattices allow a 
good control of the bandgap (in fact better than obtained in the growth of alloys) which 
can be a variable from zero to about 1.6 eV. We have also done some rough modeling 
of high mobility devices which are based on doping in the barrier to furnish electrons 
to the quantum wells. Some results are given below. Representative references are 
listed. 

There are several comments one could make about the utility of superlattices as IR 
detectors. In some ways superlattices are an advantage. For one, the energy gap 
versus SL well layer thickness d w can be better controlled in SL’s than E g versus x in 
bulk crystals - at least for E g ’s of interest. This is particularly true because the cut off 
wavelength for E g ’s of interest is a much more rapidly varying function of x for bulk 
crystals than of the layer thickness for superlattices. 

There exists a neat trick (Stormer and Dingle) in SL’s that can be used to increase the 
number of electrons without reducing the low temperature mobility due to ionized 
impurity scattering. This is accomplished by adding donors to the barrier region 
allowing the electrons to fall into the “wells”. In calculations dealing with electron 
mobility in superlattices, it has been found that the temperature dependence of the 
mobility is greatly affected by interface roughness which leads to complications in 
making accurate calculations. 

There are other advantages to the HgTe-CdTe superlattices as given by Smith, McGill 
and Schulman. The superlattice tunneling currents are reduced compared to bulk and 
diffusion currents in photovoltaic devices made from SL’s are reduced as is band to 
band tunneling currents across the junction. 

HgTe-ZnTe SL’s have also been considered but because of a 6.5% lattice mismatch 
these become strained layer superlattices. The straining can cause additional effects 
which may be tailored to advantage. We have also done calculations which model 
crudely the way mobility can be enhanced by using layered structures consisting of 
alternate layers of “quantum wells” and barriers. 
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The electron wave function in superlattices is approximately (1); 


A 

- 

r 


i k x r A -* A 

2 C e V cf (r) Xmj z ) ifr is in an A layer, 


and 


B ^ ik , r B B — 

4> - = e ^ciT ^ if r is in a B layer, 

(r) c 


where is the direction perpendicular to the z direction (perpendicular to the layers), 


AorB 


V js the Bloch wave function for the A o r B material, and Xm (z) is the envelope 
wave function which is given by the Schroedinger equation: 


h' 


+ V r (z) 


2m (z) dz 


AorB * „ 

Xm (Z) € m X m * z » rB 


I 


in a reasonable approximation. 

Here, m* is the electron effective mass of the A or B material, V c (z) is the energy level 

of the bottom of the conduction bands and e m is the energy eigenvalue which 

represents the confinement energy of the carriers. The boundary conditions at the 

interfaces between A and B materials are that X (z) and [1/m* (z) ] [3X m (z)/3z] 

m 

be continuous. The energy levels in the conduction band can be approximated by the 
Kane model for describing the electrons and holes in the A and B materials. Further 
details are given in Bastard (2). 


In general the question of energy bands in superlattices is very complicated. For one 
situation, once can picture the parabolic quarter wells as splitting the original A or B 
bands into mini bands. One of the most controversial items in band structure 
determinations is the value of the valence band offset. There is also band bending 
which can occur when charge is transferred between barriers and wells. 
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In the following graphs we have done a calculation which models very crudely the way 
that mobility can be enhanced by using layered structures consisting of alternate 
layers of quantum wells and barriers. If one adds donors in the barriers the resulting 
electrons will drift to the quantum well resulting in more electrons without more 
scattering centers. The actual calculation we have done keeps the number of 
electrons constant and reduces the number of scattering centers. Thus the low 
temperature mobility (parallel to the layers) would be increased. In the following 
graphs sit = 10 means we have electrons at low temperatures with only N^/IO 

ionized impurities/volume. Further increase of sit to 100 produces minimal low 
temperature gain in the mobility. There is no change in the effect for further increase 
in sit which seems to mean that the effect “saturates” for high electron concentration at 
low temperature. This would mean that ii impurity scattering has become negligible 
due to a decrease in ii scattering centers. 


mobility (cm2/V.s) 
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A superlattice model for MCT 
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H. MMT and Parameter Tables 

A preprint of a paper is enclosed which has the parameter tables for MMT and 
also the results of our mobility calculations. 
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ABSTRACT 

Electron mobility in the n-type diluted magnetic semiconductor 
Hgi- x Mn x Te has been calculated in the absence of a magnetic field. 

The results were compared to those of Hgi- x Cd x Te for the same range of 
energy gap. The calculations are ba sed on solving the Boltzman equation 
using a variational principles technique. The scattering processes that 
are included in the calculations are the scattering of the electrons with 
the ionizied impurities, holes, alloy disorder, acoustic phonons and 
longitudinal optical phonons. At room temperatures the mobilities are 
nearly the same, while at temperatures of order of liquid nitrogen they 
are about 30% less for MMT as compared to MCT with the same 
concentration of defects. Hgi- x Mn x Te represents a good material for 
infrared detectors. One of the advantages it has over Hgi- x Cd x Te is that 
the band gap changes faster with Mn in Hgi- x Mn x Te than it does with Cd in 
Hgi- x Cd x Te which leads to less scattering due to alloy disorder. 


) 
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I. Introduction 

Although Hgi- x Cd x Te (MCT) has many problems due to lattice, 
surface, and interface instabilities^) it is still a dominant material for 
use in infrared detectors. Because of the disadvantages of MCT other 
materials are being studied as better alternatives. One may start the 
search for new materials by noting that Cd, Zn, Mg and Mn open the band 
gap of the semimetals HgTe and HgSe( 2 ) and thus produce potential 
infrared detector materials- Hgi- x Zn x Te may be superior for use as an 
infrared detector over Hgi. x Cd x Te because it is chemically more stable 
and at the same time shares the same electronic properties^- 7 ). Hgi- 
x Cd x Se typically shows n-type conductance and does not easily show type 
conversion^). Mg has a stability problem with the Hg-Te bond when the 
required amount for the 10|im range is introduced^). Hgi- x Mn x Te (MMT) is 
not a true 1 1- VI alloy but it acts as one. The 3d orbitals of Mn are exactly 
filled with one electron per orbital and the five spins in these orbitals are 
parallel, by Hund's rule. It would require about 6-7 eV to add an addition 
electron with opposite spin to the Mn atom. As a result the 3d 5 shell acts 
as if it is complete shell and the Mn atom behaves as a group II 
element(io). 

MMT being a diluted magnetic semiconductor, represents a good 
material for infrared detectors for the following reasons : 

1) Cdi- x Mn x Te provides a suitable substrate for epitaxial growth due to 

its lattice constants 1 ). 

2) The band gap changes approximately twice as fast with Mn in MMT as 
with Cd in MCT. A smaller amount of Mn is needed to produce the same 
required energy gap which leads to less scattering due to alloy 
disorder^ 2 ). For detection in the 8-14 A* 777 window 20% CdTe in MCT 

is needed while only about 11% MnTe in MMT is required. 

3) Using a magnetic field, at very low temperatures, the energy gap can 

be fine tuned due to the presence of the exchange interaction03). 

4) Lattice, surface and interface instabilitiesO 4 ) in MMT are less when 

compared to MCT due to the bond stability of Hg-Te. 

5) Segregation(i5.i6) of CdTe in MCT is more of a problem than segregation 
of MnTe in MMT. Thus MMT is more uniform. 
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Diluted magnetic semiconductors (DMS) are characterized by the 
presence of substitutional magnetic ions which causes spin-spin exchange 
interactions between the localized magnetic moments and the band 
electrons. This interaction changes the band structure and the impurity 
states in the presence of a magnetic field. The significance of this 
interaction has been explained in several review articlesO 7 * 21 ). 

Single crystals of MMT have been grown using different techniques: 
the Bridgman method(22-24) ( the travelling-solvent method(25), the liquid- 
phase epitaxy method (LPE(26)), the two phase liquid method^) and the 
isothermal vapor-phase epitaxy method (VPE( 28 )). The as grown 
crystals(29) of MMT and MCT are often p-type due to Hg vacancies. By 
annealing(30) the samples for a certain time in Hg vapor at low 
temperature they can be converted to n type. Delves(3i) in 1965 started 
the work on semimagnetic semiconducting single crystals which Holm(23) 
and Furdyna and Leibler(32) et al. continued. Much of this research as well 
as work of more recent vintage, has been reviewed by Furdyna(33,34). a 
thorough review of MMT with more than 150 references has been given 
by Rogalski(2). 

Because MMT is paramagnetic, has the ZnS structure for all 
compositions at least up to x =0.35, and because of many other 
similarities, Johnson(29) et al. assumed that in the absence of an applied 
magnetic field the physical properties of MMT for x <0.35 are similar to 
those of MCT. This assumption was made based on the work of Furdyna(30) 

5 

who mentioned in his paper that mobilities ofA*“ '0 cm2/V.s or more 
were observed at 77K. However, there have not yet been electron mobility 
calculations for MMT which could be compared to MCT. The purpose of this 
paper is to collect the required band and material parameters and do the 
mobility calculations. Most of the parameters are obtained from 
experiment. They will be combined with a calculational model already 
used for MCT and MZT, transport properties of MMT in the absence of a 
magnetic field will be obtained, and the results will be compared with 
those of MCT for the same energy gap. 

Since the band structure of MMT in the absence of a magnetic field 
does not show any effects caused by the magnetic ions it can be described 
by Kane's theory(35,36). Under this condition the procedure of calculating 
the electron mobility in n-type MMT and the concentration of the free 
carriers would be the same as in our previous calculations of MZT(3,4,37). 
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An outline of the model used is given in section II, values of the band 
and material parameters extracted from experiments are given in section 
III, calculations and results are given in section IV and finally conclusions 
are given in section V. 

II Theory: 

As was mentioned before the wave functions of the conduction 
electrons and the dispersion relation are described by the Kane three band 
model ( 3 ) in the absence of a magnetic field. The non-parabolicity of the 
energy band has been considered. Due to the high ionicity of the material 
and the importance of the optical phonon interaction, the relaxation time 
approximation could not be used and a numerical technique has to be used 
in solving the Boltzman equation. In setting up the Boltzmann equation, 
the interaction of the electrons with the ionized impurities, heavy holes, 
alloy disorder, acoustic phonons and polar optical phonons were 
considered. At low temperature the scattering is mainly due to ionized 
impurities. The charged impurity atoms do not scatter high energy 
electrons and this scattering process decreases with increasing 
temperature. Optical phonon scattering is predominant at high 
temperature and the Frolich perturbation potential is used to describe the 
interaction. For acoustic phonons, the deformation potential which arises 
from the changes in the energy band gap as a result of the changes in the 
spacing of the lattice atoms is used. In mixed alloys, like MMT, Hg and Mn 
are randomly distributed in the crystal. Since a virtual crystal is 
assumed for band structure purposes, the disorder in the crystal is taken 
care of by the scattering of electrons from the alloy disorder. The details 
of the calculations of each scattering mechanism, the model used for the 
donor levels, the neutrality equation and the variational principles 
technique used in solving the Boltzmann equation are given in (3, 37, 38- 
42). Screening effects are included in the interaction with ionized 
impurities, holes and optical phonons. 

Since MMT contains the magnetic ions Mn++, the total electronic 
Hamiltonian(30) is given by 

H = H+H m =tl + 'ZJ(r-R l )S r < 7 , ( 1 ) 

", 

where j S the Hamiltonian term for a nonmagnetic semiconductor, H gx i s 
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the term which takes care of the exchange interaction of the band 

electrons with the localized magnetic moments, $i and <j are the spin 

operators of Mn++ and of the band electron, respectively, J is the electron- 

ion exchange coupling constant, and r and ^ are the band electron and 

Mn++ coordinates, respectively. The summation is taken over all the 

lattice sites occupied by Mn++ ions. Furdyna(3°) pointed out that since the 

electronic wave function is very extended, the electron encounters many 

Mn++ ions at any time and Si could be replaced by the thermal average 
— ^ 

< 5 > of the Mn++ spin of all ions. For a magnetic field in the z direction 

— * A 

and in the paramagnetic phase < S > can be replaced be < S z > k . Thus, 
also using the virtual crystal approximation, 

H„= a,< S 2 > x£j(r -r), (2) 

where the sum now runs over all sites. < S z > is proportional to the 
magnetization which vanishes in the absence of an external field for a 
paramagnet. It is then clear from equation (2) that the exchange term will 
also vanish in the absence of a magnetic field, at least in the 
paramagnetic phase. MMT is paramagentic( 3 °) at the composition and 
temperature we are interested in. 

III. Parameters for electron mobility and charge carrier 
concentration calculation for Hgi- x Mn x Te : 

Mercury telluride crystallizes in the zinc blende type structure with 
a lattice constant a = 6.4614°, while manganese telluride crystallizes 
in the hexagonal NiAs type structure with lattice constants! 43 ) a = 4.15 4° 
and c = 6.71 4° . X-ray measurements and differential thermal analysis 
were performed by Delves and Lewis( 44 ) for the HgTe-MnTe system and the 
following three phases were observed in the solid solution: (1) a phase 
with the zinc blende type structure (0 < x <0.36), (2) a phase plus MnTe2 
with the pyrite structure (0.3 6 < x <0.75), and (3) MnTe with NiAs type 
structure (x > 0.75). For x <0.36, Hgi- x Mn x Te crystal growth behaviour is 
very similar to that of Hgi. x Cd x Te for comparable! 3 ) x . For x >0.36, Hgi- 
x Mn x Te has more than one crystallographic phase which makes the 
material not favorable for devices applications. 

The energy gap in MMT behaves in the same manner as in the MCT. It 
changes from a negative band gap semimetal to a positive band gap 
semiconductor as the composition x, (which gives the concentration of Mn 
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in MMT or Cd in MCT), increases from zero. A zero band gap is achieved at 
4.2K temperature^), when x=0.07 in MMT and when x=0.15 for MCT. The 
change in MMT energy gap with x is twice as fast as it is in MCT(30). The 

equation deduced from experiments for the energy gap ( “ *b ) as a 

function of those two terms has been given by Kaniewksi and Mycielski(45), 

Egix, T) =-0.253 + 3. 446x +4.9 * 10"V -2.55* 10 ~*xT (eU). 

The Kane momentum matrix element was obtained from interband 
magnetooptics experiments (46,47) and it is a function of the 
composition^) x, P(x) = (8.35 - 7.94x) x ]0~°eVcm . 

Gebicki and Nazarewicz(49) measured the infrared reflection spectra 
in the optical region from 50 to 300 cm- 1 at 90 and 295 K of MMT with 
0< x <0.69. These crystals were grown by the Bridgman method and 
their crystallographic structure was tested by X-ray diffraction to ensure 
they have the zincblende structure. They obtained the dielectric functions 
versus wave number by applying the Kramers-Kronig analysis to the 
reflection spectra. The position of the maxima of Ime ( co ) and 
lm(-Ve(£o)) indicate the long-wavelength phonon frequencies for 
transverse and longitudinal optical phonons respectively. For each value 
of x, two frequencies of the transverse optical phonons and two of the 
longitudinal optical phonons were found as is consistent with "two-mode" 
behaviour" of the phonon spectra. They calculated the long-wavelength 
optical phonon frequencies as a function of composition using the random- 
element isodisplacement (REI) model (so.si) which is applicable only for 
the composition range at which the crystal structure does not change. For 
this reason the unknown parameters of the REI model were determined 
from the optical phonon frequencies for x = 0 (HgTe) and x = 0.69. The x = 

1 (MnTe) data cannot be used because the change of the crystal structure 
causes discontinuity of the phonon spectra. When they compared their 
theoretical results for the composition range 0<x<0.69 with 
experiment they obtained good agreement. Using their results for the 
dependence of long-wavelength optical phonon frequencies on composition 
we extracted the frequency equations we needed for our model using 
linear interpolation between the points x = 0 and x = 0.4 (where the curve 
was linear). We obtained: 

u LO (HgTe ) = 134-5x t © ro ( HgTe ) = 1 10 + 27. 5x 

co LO (MnTe ) = 1 87 + 25x t o> ro (MnTe) = 1 87 - 5x 
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The Callen effective charge, e’ c , which is related to the relative 
displacements of atoms in a unit cell, determines the strength of the 
interaction of the electrons with the LO-phonon (see reference 39, page 
80). The Szigeti effective charge, e ' s , and the total transverse charge are 
related to the Callen effective charge by the relations: 


3 - [e- + ‘‘ 
' C- L 3e_ 


and 


e = — 


In reference (51), Gorska et al. made the assumption that the charges ee 
and e c in the REI theoretical model are the Szigeti effective charges 


e B ~ e AB' 


e c=e AC , 


W. Gbicki et al.! 48 ), in their calculation of optical phonons in MMT 
took the effective charge of all ions to be equal to that of HgTe! 42 ) 
eA/e = ee/e = e c /e = 0.6 

Using these equations and the values mentioned earlier for the dielectric 
constants for HgTe and MnTe, the transverse charges are: 

e^(HgTe) = 3.44, ^{MnTe) = 3.32. 


The parameters needed for the interaction of the electrons with the 
acoustic phonons are the longitudinal and transverse sound velocities and 
the deformation potentials! 42 ). The longitudinal and transverse acoustic 


waves has the simple relations^) , ^ = V P and r,- V p . Because 
of the lack of measurements done on MnTe and due to the fact that the Mn 
portion is very small, we assume the Mn causes a negligible change in V| 
and vt. Using the longitudinal and Transverse sound velocities^) of HgTe, 
Cn = 5.361 x 10 11 dyne/cm 2 , C44 = 2.123 dyne/cm 2 and the density ( 2 ) of 
HgTe = 8.076 (gm/cm3), the longitudinal and transverse velocities can 
easily be estimated. We obtain V| = 2.576 x 105 cm/sec and v t = 1.621 x 
105 cm/sec. 


7 ? 



For MMT a virtual crystal was assumed and the band structure was 
described by the Kane theory. The compositional disorder is considered in 
the scattering of the electrons. It is treated as a random distribution of 
square-well scattering centers with dimensions of a unit cell and with a 
depth approximately equal to the difference between the band gaps of the 
end-point compounds! 38 ). Using the Energy gap! 2 ) values at 4.2KQ of HgTe 
and MnTe when they both have the zinc blende structure we got a value of 
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3.5eV for the disorder energy (Edis)* These results as well as values for 
all other needed quantities are summarized in Table (1). 

IV. Calculations and Results: 

MCT and MMT are compared for two different values of MMT enrgy 
gap. The first value of MMT energy gap, x (MMT) = 0.0923, is chosen to give 
values of the MMT energy gap that are very close to those of MCT, with 
x(MCT) = 0.193, at low temperature (see fig. 1). This value was chosen as 
typical for energy gaps appropriate to the 10/*^? window. Other 
appropriate values of x can be used with equal facility in our mobility 
calculations. The second value of MMT energy gap, x(MMT) = 0.094, is 
chosen to give values of the MMT energy gap that agrees with those of MCT 
with x = 0.193 at intermediate temperature (see fig. 3). In both cases we 
used a donor concentration equal to Nd = 3 x IO 1 4 /cm3 and acceptor 
concentration equal to N a = 0. In each case we considered two kinds of 
donors, one kind is ionized with an energy which merges with the 
conduction band and a concentration chosen so that Nd - N a is equal to the 
concentration of electrons in MCT at very low temperature^). The second 
kind are bound donors with a concentration and ionization energy chosen to 
fit the concentration of the electrons vs temperature at low temperature 
for MCTHO). For the concentrations considered, the second donor had very 
little effect. When the electron mobility in both cases were calculated 

and compared they came out to be very close (see graphs 2 and 4) at high 
temperatures but differed by about 30% at liquid nitrogen temperatures. 
When another value of the MMT disorder energy, 1.8eV, was used the 
mobility value was increased by approximately 25% over the whole range 
of temperature. 

V. Conclusions 

Since MMT shares the same electronic properties as MCT for the 
same energy gap and since it has the advantages that were mentioned 
earlier, it is another strong candidate for use in infrared detectors. 
However, Brimrose apparently continues to be the only company which 
makes MMT detectors commercially. The possibility of varying the 
electronic properties of MMT by use of a magnetic field makes MMT a very 
intertesting material. The effects of this on the mobility await further 
study. 
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Table 1 


i 

i 

i 


Parameters for electron mobility and charge carrier concentration 
calculations for Hgi- x Mn x Te 


Parameter value reference 


Lattice constant 

a(x) = 6.461 - 0.121x A° 

(2) 


(0.0 8 < x < 0.30) 


spin-orbit splitting 

A{eV) = 1.08 

(2) 


(0.0 8 £ x <0.30) 


r 6 - r 8 energy gap 

Eg(x t T) = -0.253 + 3.446x + 4.9 x 

(4 5) 


10"V -2.55x1 0" 3 xT (eV) 
(0.08 < x £ 0.30) 


heavy hole effective mass m 'J m o~ 0-5 (2) 

momentum matrix element P(x) = (8.35 -7.94x) x AO^eVcm (48) 

coupling conduction and valence bands 

static dielectric constant e 0 = 20.5-32.6x +25.1X 2 (2) 

high frequency dielectric constant e„ = 15.2 -28.8x +28.2X 2 (2) 

HgTe LO phonon frequency ^> LO ( H 9 Te ) = 1 34 - 5x cm -i (49) 

HgTe TO phonon frequency o>To( H 9 Te ) = 1 1 0 + 27. 5x cm -i (49) 

MnTe LO phonon frequency (o LO (MnTe ) = 1 87+ 25x cm -i (49) 

MnTe TO phonon frequency (MnTe) = 1 87 - 5x cm -i (49) 

transverse effective charges HgTe ej/e = 3.44 (39,42,49,51) 

MnTe ej/e = 3.32 

deformation potentials^) E 0 = 9.5 eV longitudinal mode (53) 

E-| ( a ) = 4 eV transverse mode (42) 

E 2 = 3 eV transverse mode (54) 

longitudinal sound velocity V| ■- 2.576 x 105 cm/sec (52) 

transverse sound velocity v t = 1.621 x 105 cm/sec (52) 

Hg mass m Hg = 200.59 a.m.u. 

Mn mass m Mn = 54.938 a.m.u. 


Te mass m Te = 127.60 a.m.u 

reduced mass m Hg Te = 1-295 xlO -22 g 
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mMnTe = 0.615 X 10' 22 g 

disorder energy Edis = 3.5eV 


a) using hydrogenic approximation 


(2), (38) 


/ 
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I. References 

We have given above a few particularly relevant references for each topic. An 
extensive list of references has been computerized by us. We have listed it by topics 
in the appendix. 
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J. Availability of Program 

We have the program written in Fortran, and stored on a Sun, a Vax and a Nighthawk. 
However, because we have local networks as well as Internet, we can move the 
program around easily from computer to computer (on campus, in state or even out of 
state). Thus we could send our program or results to a properly networked terminal at 
Marshall Space Flight Center. The programs for MCT, MZT and MMT are also 
available on disks that are available as required. 
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ABSTRACT 


The electron mobility, concentration of electrons, light holes and heavy holes, 
and Fermi energy in n-type Hg-|. x Zn x Te (MZT) have been calculated and compared 
with those of Hg.|_ x Cd x Te (MCT) for the same energy gap (but different x), and the 
same donor and acceptor concentration. The results for MZT were found to be very 
close to those of MCT. The calculation relied on solving the Boltzmann transport 
equation using variational principals. The processes that were of significance 
included the scattering of electrons by: ionized impurities, holes, compositional 
disorder, acoustical phonons and optical phonons. In the process of calculating the 
results for MZT a table for the latest values of all the needed material parameters was 
obtained by a literature search. 

MZT has been considered as a substitute for MCT for possible use in infrared 
detectors due to its relative stability and hardness. Our theoretical results for MZT 
were compared to the experimental results published by Rolland et al. and Granger et 
al. and very good agreement over most of the temperature range was obtained. 

The effect of making changes on some of the relevant parameters in the 
calculation was checked. Also, adding neutral impurities to the scattering processes 
and using the most recent published band parameters of MCT (rather than the slightly 
older values) on the electron calculation did not produce significant changes. 

The speed of a very involved and long program has been increased by 
changing its code from Basic for an HP 9845 to FORTRAN for a Vax 1 1/780 and later a 
Sun workstation 386i. This program was given to us by Dr. S.L. Lehoczky and it was 
written to calculate the electron mobility, the free carrier concentrations, and the Fermi 
energy for MCT. All of these were calculated as a function of temperature, 
composition, donor, and acceptor concentration. The program was checked, corrected 
and fine tuned by including the latest values of material parameters for MCT, another 
version was made to fun for MZT and neutral defect scattering was added to the 
program. 


